APPROXIMATING AMOEBAS AND COAMOEBAS 
BY SUMS OF SQUARES 



THORSTEN THEOBALD AND TIMO DE WOLFF 

Abstract. Amoebas and coamoebas are the logarithmic images of algebraic varieties 
and the images of algebraic varieties under the arg-map, respectively. We present new 
techniques for computational problems on amoebas and coamoebas, thus establishing 
new connections between (co-)amoebas, semialgebraic and convex algebraic geometry 
and semidefinite programming. 

Our approach is based on formulating the membership problem in amoebas (respec- 
tively coamoebas) as a suitable real algebraic feasibility problem. Using the real Null- 
stellensatz, this allows to tackle the problem by sums of squares techniques and semi- 
dchnite programming. Our method yields polynomial identities as certificates of non- 
containedness of a point in an amoeba or coamoeba. As main theoretical result, we 
establish some degree bounds on the polynomial certificates. Moreover, we provide some 
actual computations of amoebas based on the sums of squares approach. 



1. Introduction 

For an ideal / C C[Z\, . . . , Z n ] an amoeba (introduced by Gel'fand, Kapranov, and 
Zelevinsky [3], see also the surveys [8] or [16]) is the image of the variety V(J) where 
each complex coordinate is mapped to (the logarithm of) its absolute value. It is often 
customary and useful to consider the logarithmic version of an amoeba 

Ar = {(log N,..., log \z n \) : zGV(7)n(C*n 

with C* := C \ {0}, and we denote the unlog amoeba by 

Uj = {(M,..., |^|) : ZEV(I)}. 

Similarly, the coamoeba Ci is defined as Cj := Arg(V(J) PI (C*) n ), where Arg denotes the 
mapping 

(zi, ...,z n )^ (arg(*i), . . . , arg(z n )) e T n := (M/2vrZ) n 

and arg denotes the argument of a complex number (see [TUl [III [121 [13]). If / is a principal 
ideal generated by a polynomial /, we shortly write Af := Am and analogously Uf, Cf. 

Studying computational questions of amoebas has been initiated in [24], where in par- 
ticular certain special classes of amoebas (e.g., two-dimensional amoebas, amoebas of 
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Grassmannians) were studied. One of the natural and fundamental computational ques- 
tions is the membership problem, which asks for a given point (Ai, . . . , A n ) whether this 
point is contained in an (unlog) amoeba respectively coamoeba. 

In [19] Purbhoo provided a characterization for the points in the complement of a 
hypersurface amoeba which can be used to numerically approximate the amoeba. His 
lopsidedness criterion provides an inequality-based certificate for non-containedness of 
a point in an amoeba, but does not provide an algebraic certificate (in the sense of 
a polynomial identity certifying the non-containedness). The certificates are given by 
iterated resultants. With this technique the amoeba can be approximated by a limit 
process. The computational efforts of computing the resultants are growing quite fast, 
and the convergence is slow. 

A different approach to tackle computational problems on amoebas is to apply suitable 
Nullstellen- or Positivstellensatze from real algebraic geometry or complex geometry. For 
some natural problems a direct approach via the Nullstellensatz (applied on a realification 
of the problem) is possible. Using a degree truncation approach, this allows to find sum- 
of-squares-based polynomial identities which certify that a certain point is located outside 
of an amoeba or coamoeba. In particular, it is well known from recent lines of research 
in computational semialgebraic geometry (see, e.g., [UEKE]) that these certificates can 
be computed via semidefinite programming (SDP). 

In this paper, we discuss theoretical foundations as well as some practical issues of 
such an approach, thus establishing new connections between amoebas, semialgebraic 
and convex algebraic geometry and semidefinite programming. Firstly, we present various 
Nullstellensatz-type formulations (Statements 13.21 and 13.71) and compare their properties 
to a recent toric Nullstellensatz of Niculescu and Putinar [9J. Using a degree trunca- 
tion approach this yields effective approximation hierarchies for the amoeba membership 
problem (Theorem 13.111) . 

The main theoretical contribution is contained in Section HI For one of our approaches, 
we can provide some degree bounds for the certificates (Corollary 14.51) . It is remark- 
able and even somewhat surprising that these degree bounds are derived from Purbhoo's 
lopsidedness criterion (which is not at all sum-of-squares-based). We also show that in 
certain cases (such as for the Grassmannian of lines) the degree bounds can be reduced 
to simpler amoebas (Theorem 14.71) . 

In Section Owe provide some actual computations on this symbolic- numerical approach. 
Besides providing some results on the membership problem itself, we will also consider 
more sophisticated versions (such as bounding the diameter of a complement component 
for certain classes). 

Finally, in Section [6] we give an outlook on further questions on the approach initiated 
in the current paper. 

2. Preliminaries 

In the following, let C[Z] = C[Z X , . . . , Z n ] denote the polynomial ring over C in n 
variables. For / = XLeA c a^ a e QZ], the Newton polytope New(/) = conv{a e Nq : 
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Figure 1. The amoeba of / = X 2 Y + XY 2 - AXY + 1. 



a G supp(/)} of / is the convex hull of the exponent vectors, where supp(/) denotes the 
support of /. 

2.1. Amoebas and coamoebas. We recall some basic statements about amoebas (see 
[21 US]). For any ideal / C C[Z], the amoeba A\ is a closed set. For a polynomial 
/ G C[Z], the complement of the hypersurface amoeba Af consists of finitely many 
convex regions, and these regions are in bijective correspondence with the different Laurent 
expansions of the rational function l/f. See Figured] for an example. 
The order v of a point w in the complement MJ 1 \ Af is given by 



where Log(z) is defined as Log(z) = (log \zi\, . . . , log \ z n \). The order mapping induces an 
injective map from the set of complement components into New(/) fl27\ The complement 
components corresponding to the vertices of New(/) do always exist [2]. 

Similarly, for / G C[Z] any connected component of the coamoeba complement T n \ Cf 
is a convex set (see Figure [2]). If Cf denotes the closure of Cf in the torus T n , then the 
number of connected components of T n \ Cf is bounded by n! volNew(/) (Nisse [131 Thm. 
5.19]), where vol denotes the volume. 

For technical reasons (as explained below) it will be often convenient to consider in 
the definition of a coamoeba also those points z G V(I) which have a zero-component. 
Namely, if a zero z of / has a zero- component Zj = then we associate this component 
to any phase. Call this modified version of a coamoeba Cj. Note that for principal ideals 
/ = (/) the difference between Cj and C\ solely may occur at points which are contained 
in the closure of Cj. The set-theoretic difference of Cj and C\ is a lower- dimensional subset 
of M. n (since in each environment of a point in C\ \ Cj we have a coamoeba point). 




Tl 



±<J<n, 
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Figure 2. The coamoeba of / = X 2 Y + XY 2 - AXY + 1 in two different 
views of T 2 , namely [— 7r,7r) 2 versus [0,27r) 2 . 

2.2. The situation at oo. It is well-known that the geometry of amoebas at infinity 
(i.e., the "tentacles") can be characterized in terms of logarithmic limit sets and tropical 
geometry, and thus amoebas form one of the building blocks of tropical geometry (for 
general background on tropical geometry we refer to [3 (HJ [21] ) . 

For (large) R > let denote the scaled version A\ R ^ := j^Ai fl S" --1 , where § n_1 
denotes the (n — l)-dimensional unit sphere. Extending this definition, the logarithmic 
limit set A^ is the set of points v G §™ _1 such that there exists a sequence vr G ,4^ 
with linifl-^oo Vr = v. For a polynomial / = J2 a c aZ a G C[Z] denote by trop/ : = 
00 Z^ 1 • • • Z^ n its tropicalization (with respect to the trivial valuation) over the 
tropical semiring (R, ©, 0) := (R U {— oo}, max, +). Then (see [7J [23]): 

Proposition 2.1. A vector w G lR n \ {0} is contained in the tropical variety of I if and 
only if the corresponding unit vector rr^nW is contained in Wf°\ Thus the tropical variety 

of I coincides with the cone over the logarithmic limit set A^. 

3. Approximations based on the real Nullstellensatz 

We study certificates of points in the complement of the amoeba based on the real 
Nullstellensatz and compare them to existing statements in the literature (such as the 
toric Nullstellensatz of Niculescu and Putinar). By imposing degree truncations this will 
then yield a hierarchy of certificates of bounded degree. 

We use the following real Nullstellensatz (see, e.g., [U ITS]): 

Proposition 3.1. For polynomials gx, g r G R[X] and I := (gi,...,g r ) C M.[X] the 
following statements are equivalent: 

• The real variety Vr(7) is empty. 
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• There exist a polynomial G G / and a sum of squares polynomial H with 

G + H+l = 0. 

Given A G (0, oo) n , the question if A is contained in the unlog amoeba Ui can be 
phrased as the real solvability of a real system of polynomial equations. For a polynomial 
/ G C[Z) = C[Z U ...,Z n ) let r, f m G R[X, Y] = R[X U . . . , X n , Y u . . . , Y n ] be its real 
and imaginary parts, i.e., 

f(Z) = f{X + iY) = f e (X,Y)+i-r(X,Y). 

We consider the ideal I' C R[X, Y] generated by the polynomials 

(3.1) {./J'- .//" : l<J<r}u{X 2 k + Y k 2 -Xl : 1 < k < n) . 

Corollary 3.2. Let I = (fi,---,f r ), and A G (0, oo) n . Either the point A is contained 
in Hi, or there exist a polynomial G G I' C M.[X, Y] and a sum of squares polynomial 
H G R[X, Y] with 

(3.2) G + H+l = 0. 

Proof. For any polynomial / G C[Z] it suffices to observe that a point z = x + iy is 
contained in V(/) if and only if (x,y) G V]R(/ re ) PI Vu(f im ), and that \z^\ = Xk if and only 
(x, y) G Vr(X 2 + Y 2 — X\). Then the statement follows from Proposition 13.11 □ 

Corollary 13.21 states that for any point X^LlAi there exists a certificate 

r r n 

j=i j=i k=i 

with polynomials Pj,p'j, qk and a sum of squares H. 

Remark 3.3. By the following fact (which is easy to check), the sum of squares condi- 
tion EH] can also be stated shortly as 

—1 is a sum of squares in the quotient ring R[X, Y]/I' . 

Fact 3.4. (Parrilo [15J.) Let I = (g u ...,g r ) C R[X] and f G R[X\. There exist 
Pi, ■ ■ ■ ,Pk £ R[X] such that 

f + ^^Pigi is a sum of squares in R[X] 

i 

if and only if f is a sum of squares in R[X]/I. 

Any two of these equivalent conditions in Fact 13.41 is a certificate for the nonnegativity 
of / on the variety /. 

Before stating a coamoeba version, we note the following normalization properties. 
Whenever it is needed for amoebas, we can assume that A is the all- 1- vector 1. Similarly, 
for amoebas we can assume that all components of the point under investigation has 
phase 0. 
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Lemma 3.5. Let I = (fi, . . . , f r ). 

(1) A point (Ai, . . . , A n ) G (0, oo) n is contained inUi if and only iflE M(gi,...,g r ) , where 

g j (Z 1 ,...,Z n ) := f j (X 1 Z 1 ,...,X n Z n ), l<j<r. 

(2) A point (zi, . . . , z n ) is contained in V(I) with argz,,- = fij if and only if the (non- 
negative) real vector y with yj := Zje~ tflj is contained in V(gi, . . . ,g r ) where 

9j {Z u ...,Z n ):= f 5 {Z x e^\. . . , Z n e^) , 1 < j < r . 

Proof. A point (z\, . . . , z n ) is contained in V(J) with \zj\ = Xj if and only if the vector y 
defined by yj := Zj/Xj is contained in V(g±, . . . ,g r ) with \yA = 1. The second statement 
follows similarly. □ 

Theorem 3.6. Let I = (f 1 , . . . , f r ) . The point (0, . . . , 0) is contained in the complement 
of the coamoeba C\ if and only if there exists a polynomial identity 

r r n 

(3.3) c * ■ ^ X2 ' F ) r ° + E < ■ ^ X2 > + zZ d r Y 3 + H + l = 

i=i i=i j=i 

with polynomials Ci, dj G M[X, Y] and a sum of squares H. Here, fi(X 2 , Y) abbreviates 
/ i (X 1 , . . . , X n , Yx, . . . , Y n ) . 

Proof. Note that the statement G C\ is equivalent to {z = x + iy G C n : z G 
V(I) and Xi > 0, yi = 0, 1 < i < n} ^ 0. Moreover, observe that the condition x,i > 
can be replaced by considering Xf in the arguments of fx, . . . , f r . Hence, by Proposi- 
tion [3J] the statement ^ C\ is equivalent to the existence of a polynomial identity of the 
form ([33]). □ 

Observe that in the proof the use of C'j (rather than Cj) allowed to use the basic 
Nullstellensatz (rather than a Positivstellensatz, which would have introduced several 
sum of squares polynomials). 

The following variant of the Nullstellensatz approach will allow to obtain degree bounds 
(see Section H]). For vectors a(l), . . . , a(d) G Nq and coefficients bi,...,bd G C* let 
/ = EtA- ■ Z a{l) eC[Z}. For any given values of Ai, . . . , A n set 



X a(j) = X «U)l . . . X *(j) n ^ l<j< d . 



If the rank of the matrix with columns a(l), . . . , a(d) is n (i.e., the vectors a(l), . . . , a(d) 
span R n ) then the A-values can be reconstructed uniquely from the /i-values. We come 
up with the following variant of a Nullstellensatz. 

Let rrij be the monomial rrij = Z a ^ = Z"^ 1 ■ ■ ■ Zn^ n ■ We consider the ideal /* C 
M[X, Y] generated by the polynomials 

(3.4) {f» , ff : 1 < j < r} U {(mf) 2 + (mf ) 2 - /if : 1 < j < d} . 
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Corollary 3.7. Let f = Y%=\h ' za{i) Wlth Qt(l), . . . , a(d) spanning R n and I = (/). 
Further let A G (0,oo) n . Either a point (Ai,...,A n ) is contained in Uj, or there exist 
polynomials G G I* C ~K[X, Y] and a sum of squares polynomial H G M.[X, Y] with 

(3.5) G + H+l = 0. 

For hypersurface amoebas of real polynomials, the membership problem relates to the 
following statement of Niculescu and Putinar [9]. Let p = p(X,Y) G M[Xx, . . . , X n , 
Yi, . . . , Y n ] be a real polynomial. Then p can be written as a complex polynomial p(X, Y) = 

P(Z, Z) with P G C[Zi ...,Z n ,Zi,..., Z n ] and P(Z,Z) = P(Z, Z). Note that there exists 
a polynomial Q G C[Zi, . . . , Z n ] with 

p(x,?/) 2 = \P(z,z)\ 2 = \Q{z)\ 2 for zeT n , 

where T := {z E C : |z| = 1}. 

The following statement can be obtained by applying the Nullstellensatz on the set {z = 
(x,y) : |5 , (^)| 2 = 1, \zi\ 2 = 1, . . . , \z n \ 2 = 1}, then applying Putinar's Theorem [20] on the 
multiplier polynomial of |g(2')| 2 (see [9]). 

Proposition 3.8. Let q G C[Zi, . . . , Z n ]. Then q(z) ^ for all z G T n if and only if 
there are complex polynomials Pi, ■ ■ ■ ,Pk, ri, . . . , r t G C[Zi, . . . , Z n ] 

(3.6) l + \pi(z)\ 2 + --' + \pk(z)\ 2 = \q(z)\ 2 (\r 1 (z)\ 2 + --- + \r l (z)\ 2 ), zeT n . 

Note that the statement is not an identity of polynomials, but an identity for all z in 
the n-torus T n . 

While Proposition 13.81 provides a nice structural result, due to the following reasons 
we prefer Corollary 13.21 for actual computations. In representation (13. 6p . two sums of 
squares polynomials (rather than just one as in (13.21) ) are needed in the representation, 
and the degree is increased (by the squaring process). Moreover, the theorem is not really 
a representation theorem (in terms of an identity of polynomials), but an identity over 
T n ; therefore in order to express this computationally, the polynomials hidden in this 
equivalence (i.e., the polynomials 1 — \Zi\ 2 , . . . , 1 — \Z n \ 2 ) have to be additionally used. 

SOS-based approximations. By putting degree truncations on the certificates, we can 
transform the theoretic statements into effective algorithmic procedures for constructing 
certificates. The idea of degree truncations in polynomial identities follows the same prin- 
ciples of the degree truncations with various types of Nullstellen- and Positivstellensatze 
in [U El H3]. It is instructive to have a look at two simple examples first. 

Example 3.9. Let / be the polynomial / = Z + zq with a complex constant z = x^ + iyo. 
The ideal / of interest is defined by 

ht := f re = X + x , 
h 2 := f im = Y + y , 
h 3 := X 2 + F 2 -A 2 . 

For values of A > which correspond to points outside the amoeba (i.e., A 2 ^ x^ + y^), we 
have Vc(-0 = and thus the Grobner basis G of {hi, h 2 , h 3 ) is G = {!}. The corresponding 
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multiplier polynomials pi to represent 1 as a linear combination XliPi^i are 
-X + xp -Y + y 1 

Hence, in particular, —1 can be written as a sum of squares in the quotient ring K[X]/J. 
The necessary degree with regard to equation (13.21) is just 2. 

For A 2 = a 2 +b 2 , the Grobner basis (w.r.t. a lexicographic variable ordering with X y Y) 

is 

X + a, Y + b. 

The point (—a, —b) is contained in Vr(I); thus in this case there does not exist a Nullstellen- 
type certificate. 

Example 3.10. Consider the polynomial / = Z\ + Z 2 + 5 with Zj = Xj + iYj. The ideal 
/ of interest is defined by 

hx := X 1 + X 2 + 5, 
h 2 := Y, + Y 2 , 
h := Xl + Y?-\l, 
h 4 := X% + Y 2 2 — \\ . 

Consider Ai = 2, Ai = 3. Using a lexicographic ordering with X\ y X 2 y Y\ y Y 2 , a 
Grobner basis is 

Y 2 \ Y, + Y 2 , X 2 + 3, Xi + 2. 

The standard monomials are 1 and Y 2 . It must be that —1 is a sum of squares in the 
quotient ring. 

Using the degree truncation approach for sums of squares we can define the set Ck as 
the set of points in M™ such that there exists a certificate, say, in the standard approach, 
of degree at most 2k. And similarly for coamoebas, where we denote the sequence by D^. 
These sequences are the basis of the effective implementation (see Section [5]). 

Theorem 3.11. Let I = (fx,...,f r ) and k := maxj |~deg fi /2] . The sequence (Ck)k>k 
converges pointwise to the complement of the unlog amoeba Ui, and it is monotone in- 
creasing in the set-theoretic sense, i.e. Ck C Ck+i for k > k . 

Similarly, the sequence (-Dfc)fc>fc converges pointwise to the complement of the coamoeba 
Ci, and it is monotone increasing in the set-theoretic sense, i.e. Dk C D^+i for k > ko. 

Proof. For any given point z in the complement of the amoeba there exists a certificate of 
minimal degree, say d. For k < \d/2] the point z is not contained in Ck and for k > \d/2] 
the point z is contained in Ck- In particular, the relaxation process is monotone increasing. 
And analogously for coamoebas. □ 

It is well-known (and at the heart of current developments in optimization of polyno- 
mial functions, see [H Q3] or the survey [S]) that SOS conditions of bounded degree can 
be phrased as semidefinite programs. Finding an (optimal) positive semidefinite matrix 
within an affine linear variety is known as semidefinite programming (see e.g. [25] for a 
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comprehensive treatment). Semidefinite programs can be solved efficiently both in theory 
and in practice. 

Precisely, any sum-of-squares polynomial H can be expressed as MQM T , where Q 
is a symmetric positive semidefinite matrix (abbreviated Q >z 0) and M is a vector of 
monomials. 

Similarly, by the degree restriction the linear combination in (13. 2p or (13. 5 j) can be 
integrated into the semidefinite formulation by a comparison of coefficients. 

4. Special certificates 

For a certain class of amoebas, we can provide some explicit classes of Nullstellensatz- 
type certificates. As a first warmup-example, we illustrate some ideas for constructing 
special certificates systematically for linear amoebas in the standard approach. Then we 
show how to construct special certificates for the monomial-based approach. 

In this section we concentrate on the case of hypersurface amoebas. 

Linear amoebas in the standard approach. Let / = aZi + bZ 2 + c be a general linear 
polynomial in two variables with real coefficients a,b,c G R. We consider certificates of 
the form (13.21) based on the third binomial formula, where we use the sums of squares 
(Xi — X 2 ) 2 and (Yi — Y 2 ) 2 . For simplicity assume a, b > 0. The representation 

(aXt + bX 2 - c)(aXi + bX 2 + c) + (aY 1 + bY 2 )(aY 1 + bY 2 ) 

-{a 2 + ab)(Xl + Y 2 - A 2 ) - (b 2 + ab)(X 2 + Y 2 2 - A 2 ) + ab(X 1 - X 2 f + ab{Y 1 - Y 2 ) 2 
simplifies to 

(4.1) (a 2 + ab)\\ + (b 2 + ab)\ 2 2 - c 2 . 

Assume that the point (Ai,A2) is not contained in the unlog amoeba Uf. In order to 
obtain the desired polynomial identity (13. 2p certifying containedness in the complement 
of Uf, we need this term to be smaller then zero. Then, by scaling, we can bring this 
to —1, so that by adding of +1 (i.e., replacing '0' in the formula by '1'), we obtain the 
polynomial identity (13.21) . 

Example 4.1. Let a = 1, b = 2, c = 5. The curve (in Ai,A2) given by (a 2 + ab)\\ + 
(b 2 + ab)X 2 — c 2 has a logarithmic image that is shown in Figure [31 By projective sym- 
metry, analogous special certificates can be obtained within the two other complement 
components. 

Since (by homogenizing the polynomial) there is a symmetry, we obtain similarly an 
approximation of the other two complement components. Hence we have: 

Lemma 4.2. For the points in the regions defined by the outer regions of the curves (14.11) 
there exist certificates of degree at most 2. 
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Figure 3. The boundary of an amoeba of a linear polynomial (blue) and 
the boundary of the outer region for which the special certificates of degree 2 
exist (red). 

The monomial-based approach. For the monomial-based approach based on Corol- 
lary I4.5I we can provide special certificates for a much more general class. Our point of 
departure is Purbhoo's lopsidedness criterion [19] which guarantees that a point belongs 
to the complement of an amoeba Af. In particular, we can provide degree bounds for 
these certificates. 

In the following let a(l), . . .,a(d) G N£ span W 1 and / = Yf j= i b j Za{i) e C[Z] with 

monomials rrij := Z a ^ = Z° ' • • • Z% . For any point v G M n define f{v} as the 
following sequence of numbers in M>o, 

f{v} := (\m 1 (Log- 1 (v))\,...,\m d (Log-\v))\) . 

A list of positive real numbers is called lopsided if one of the numbers is greater than the 
sum of all the others. We call a point lopsided, if the sequence f{v} is lopsided. 

Furthermore set 

CA(f) := {v G W 1 : f{v} is not lopsided}. 

It is easy to see that Af C CA(f) with Af ^ CA(f) in general. In the following way the 
amoeba Af can be approximated based on the lopsidedness concept. For r > 1 let 

r— 1 r— 1 

f r (Z) := J] J] / (e 2 ^ r Z 1 , e 2 ^' r Z n ) 

fc 1= k„=0 

= Res (Res (. . . Res(/(f/ 1 Z 1 , . . . , U n Z n ),U[ - 1), . . . , TT n _ x - l) , U r n - l) , 

where Res denotes the resultant of two (univariate) polynomials. Then the following 
theorem holds. 

Theorem 4.3. (Purbhoo [HI Theorem 1]) 

(a) For r — > oo the family CA(f r ) converges uniformly to Af. There exists an integer 
N such that to compute Af within e > 0, it suffices to compute CA{f r ) for any 
r > N. Moreover, N depends only on e and the Newton polytope (or degree) of f 
and can be computed explicitly from these data. 
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(b) For an ideal I C C[Z] a point v G W 1 is in the amoeba Ai if and only if g{v} is 
not lopsided for every g G /. 



Based on Theorem 14.31 one can devise a converging sequence of approximations for 
the amoeba. Note, however, that the lopsidedness criterion is not a Nullstellensatz in a 
strict sense since it does not provide a polynomial identity certifying membership in the 
complement of the amoeba. 

The aim of this section is to figure out how our SOS approximation is related to the 
lopsidedness and transform the lopsidedness certificate into a certificate for the Nullstel- 
lensatze presented in Section [31 

By Lemma [3.51 we can assume that the point w under investigation is the all- 1- vector 
1 . In this situation lopsidedness means that there is an index j G { 1 , . . . , d} with \bj\ > 
If th e lopsidedness condition is satisfied in w, then the following statement 
provides a certificate of the form G + H + 1 and bounded degree. 

Corresponding to the definition of I* in (13 .4p . let the polynomials si, . . . , s d +2 be defined 

by 

Sl = (^■(^ W ) rC ) 2 +(ff-(^ (t) ) im ) 2 -l> l<i<d, 
and s d+1 = f IC , s d+2 = / im . 

Theorem 4.4. // the point A = 1 is contained in the complement ofUf with f{0} being 
lopsided and dominating element |mi(l)|, then there exists a certificate of total degree 
2 ■ tdeg(/) which is given by 

d+2 



(4.2) J2 s i9i + H+1 = 

where 

g x = \h\ 2 , g { = -\bi\ -jTlhl, 2<i<d, 

k=2 

d \ re i ^ \ 

-h-Z^+J^bi-Z^A , g d+2 = (-h-Z^ + J^bi-Z^A 

i=2 J \ i=2/ 



i=l 



d 



9d+l 



2<i<j<d \\ *\ I 3\ / 

+w .i % i.^.(z«»r-| ;i .(z-B ) ^ . 



2 
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Proof. By the binomial theorem (a + b) ■ (a — b) = a 2 — b 2 with o,6gR, substituting the 
polynomials Si and gj into s d +igd+i + s d+2 gd+2 yields 

2 





' \ i=2 , 

Adding g^i and the SOS term H yields 

-N^ + (|:(S-(- B r) 2 + (g 

Hence, the expression Ylt=i s i9i + H in (14.21) in total results in 

-N 2 +(l> 

\i=2 

Since all \bi\ > and by our assumption of lopsidedness with dominating term |mi(l)|, 
this is the certificate we wanted to obtain. By rescaling, we can bring the constant to 
-1. □ 

We say that there exists a certificate for a point w in the complement of the (log) 
amoeba Af if there exists a certificate for the point 1 in the complement of the amoeba 
U g in the sense of Theorem 14.41 where g is defined as in Lemma T3.5I and \:—\ log" 1 (wj)|. 

Corollary 4.5. Let r G N. 

(1) For any w £ R n \ CA(f r ) C M. n \ Af there exists a certificate of degree at most 
2 • r n ■ deg(/) which can be computed explicitly. 

(2) The certificate determines the order of the complement component to which w 
belongs. 

Proof. By definition of g, we have w G Af if and only if 1 6 W 9 . Further 1 belongs to 
CA(g r ) if and only if g r {0} is not lopsided. Applying Theorem !4.4l on the function g r yields 
a certificate for w in the log amoeba Af. Since we have tdeg(^ r ) = tdeg(g)-r n = tdeg(/)-r n 
due to the definition of g r and g the result follows. 

For the second statement, note that passing over from / to g does not change the 
order v of any point in the complement of the amoebas. Now it suffices to show that the 
dominating term (which occurs in a distinguished way in the certificate) determines the 
order of the complement component. The latter statement follows from Purbhoo's result 
that if w CA(f r ) and the order of the complement component w belongs to is a(i) then 
the dominant term in f r has the exponent r n ■ a(i) (see [191 Proposition 4.1]). □ 

Corollary 4.6. For linear hyperplane amoebas in MJ 1 , any point in the complement of the 
amoeba has a certificate whose sum of squares is a sum of squares of affine functions. 

Proof. By the explicit characterization of linear hyperplane amoebas in [21 Corollary 4.3], 
any point in the complement is lopsided. Hence, the statement follows from Theorem 14.41 

□ 
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Simplified expressions. From a slightly more general point of view, the monomial- 
based certificates can be seen as a special case of the following construction. Whenever 
the defining polynomials of a variety origin from simpler polynomials with algebraically 
independent monomials, then the approximation of the amoeba can be simplified. 

For an ideal I let V := V(I) C (C*) n be its subvariety in (C*) n . Let 71,..., 7* 

be k monomials in n variables, say, 7^ = Y a ^ = Y^^ 1 Y 2 a ^ 2 ■ ■ - Yn^ n , where a(i) = 
(a(i) 1 , . . . ,a(i) n ) G Z n . They define a homomorphism 7 of algebraic groups from (C*) n 
to (C*) fc . For any subvariety W of (C*) fc , the inverse image / y~ 1 (W) is a subvariety of 
(C*) n . Note that the map 7 is onto if and only if the vectors a(l), . . . , a(k) are linearly 
independent (see [2"lj Lemma 4.1]). 

Let J be an ideal with V(J) = 7 _1 (V). If the map 7 is onto, then computing the 
amoeba of J can be reduced to the computation of the amoeba of /. Let 7' denote the 
restriction of 7 to the multiplicative subgroup (0, oo) n . Then the following diagram is a 
commutative diagram of multiplicative abelian groups: 

(C*) n (C*) fe 
(0,oo) n ^ (0,oo) fc 

where the vertical maps are taking coordinate-wise absolute value. For vectors p = 
(pi, . . . ,p n ) in (C*) n we write \p\ = (\pi\, ■ ■ ■ , \p n \) G (0, oo) n , and similarly for vectors of 
length k. Further, for V C (C*) n let \V\ := {|p| : p G V}. If the map 7 is onto then 
|7"WI =7 , ~ 1 (M) (see H). 

Theorem 4.7. // a point outside of an unlog amoeba Ui has a certificate of total degree 
d then a point outside of the unlog Uj has a certificate of degree d ■ D, where D is the 
maximal total degree of the monomials 71, . . . , 7^. 

In particular, this statement applies to the certificates from Statements 13.21 and 13.71 

Proof. Let p be a point outside of the unlog amoeba of V C (C*) n which has a certificate 
of total degree d. By Corollary (13. 2p . the certificate consists of a polynomial G(X,Y) 
in the real ideal I' C M[X, Y] from (13. ip and by real sums of squares of polynomials in 
Y]. For the polynomials in the ideals, we observe that the realification process carries 
over to the substitution process. W.l.o.g. we can assume that 7« is a product of just two 
factors. Then, with Z = Z\+ Z = PQ we have Z\ + 1Z2 = (P\ + iP2){Q\ + and 
use the real substitutions Z\ = P\Q\ — P2Q2, Z 2 = P1Q2 + PiQ\- And in the same way 
the real sum of squares remain real sums of squares (of the polynomials in P i} Qj) after 
substituting. □ 

Example 4.8. Let Gi^ denote the Grassmannian of lines in 3-space is the variety in P^, 
defined by 

-P01-P23 ~ -P02-P13 + -P03-P12 = , 

which we consider as a subvariety of (C*) 6 . The three terms in this quadratic equation 
involve distinct variables and hence correspond to linearly independent exponent vectors. 
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Figure 4. SOS certificates of linear amoebas restricted to degree 2. The 
red (dark) points represent infeasible SDPs, and in the green (light) points 
numerical instabilities were reported in the computations. 

Note that Gi )3 equals 7" 1 (V) where 

7 : (C*) 6 ->• (C*) 3 , (p l,P02,P03,Pl2,Pl3,P23) ^ (P0lP23, P02Pl3, P03P12) 

and V denotes the plane in 3-space defined by the linear equation X — Y + Z = 0. Since 
by Theorem 14.61 any point in the complement has a certificate of degree 2, Theorem 14.71 
implies that every point in the complement of the Grassmannian amoeba has a certificate 
of degree 4. 

5. Computing the relaxations via semidefinite programming 

We close the paper by providing some computational results in order to confirm the 
validity of our approach. The subsequent computations have been performed on top 
of SOSTools [17J which is a Matlab package for computing sums of squares based 
computations. The SDP package underlying SOSTools is SeDuMi [22]. 

Example 5.1. For the test case of a linear polynomial / := Z\ + 2Z 2 + 3, the boundary 
contour of the amoeba Af can be explicitly described, and it is given by the curves 

exp(^i) = 2 • exp(z 2 ) + 3 , 
2-exp(z 2 ) = exp(^ 1 ) + 3, 

3 = exp(zi) + 2 ■ exp(z 2 ) , 

see [2]. We compute the amoeba of / with our SDP via SOSTools on a grid of size 
250 x 250 lattice points in the area [— 3,4] 2 . In the SDP we restrict to polynomials of 
degree 2. By Theorem 14.41 the approximation is exact in that case (up to numerical 
issues). Figure H] visualizes the SDP-based computation of the SOS certificates. In the 
figure, at the white outer regions, certificates are found. At the red points the SDP is 
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infeasible; at the green points, no feasibility is proven, but numerical instabilities are 
reported by the SDP solver. The issue of numerical stability of our SOS-based amoeba 
computations and of general SOS computations is an important issue in convex algebraic 
geometry which deserves further study. 

Example 5.2. As in Figured] we consider the class of polynomials / := Z\Z 2 + Z\Z\ + 
c ■ Z\Z 2 + 1 with some constant c G R. We use the monomial-based approach from 
Corollary 13.71 In order to compute whether a given point (fi\, fi 2 ) G R>o is contained in 
the unlog amoeba Uf, we have to consider the polynomials 

h x = (X 2 X 2 - 2X X Y X Y 2 - Y 2 ■ X 2 ) + (JJ^Xf - X X Y 2 - 2Y 1 X 2 Y 2 ) 

+c-(X 1 X 2 -Y 1 Y 2 ) + l, 
h 2 = (X 2 Y 2 + 2X X Y X X 2 - Y*Y 2 ) + (2X 1 X 2 Y 2 + Y X X\ - Y{Y 2 ) + c • {X X Y 2 + X 2 Y X ) , 
h 3 = (X 2 X 2 -2X 1 Y 1 Y 2 -Y 2 X 2 ) 2 + (X 2 Y 2 + 2X 1 Y 1 X 2 -Y 2 Y 2 ) 2 -( f i 2 1 - f i 2 ) 2 , 
h 4 = (X 1 X 2 -X 1 Y 2 -2Y 1 X 2 Y 2 ) 2 + (2X 1 X 2 Y 2 + Y 1 X 2 -Y 1 Y 2 2 ) 2 -(fi 1 - f i 2 2 ) 2 , 
h = (X 1 X 2 -Y 1 Y 2 ) 2 + (X 1 Y 2 + X 2 Y 1 ) 2 -( f i 1 -fi 2 ) 2 . 

For the case c = 2 and c = —4 we investigate 160 x 160 points in the area [—4, 4] x [—5, 3] 
and restrict the polynomials multiplied with the constraints to degree 3. The resulting 
amoeba Af is depicted in Figure [5j At the white points the SDP is feasible and thus these 
points belong to the complement component. At the orange points the SDP is recognized 
as feasible with numerical issues (within a pre-defined range). At the black points the 
SDP was infeasible without and at the turquoise points with numerical issues reported. 
At the red points the program stopped due to exceeding numerical problems. The union 
of the black, the turquoise and part of the orange points provides the (degree bounded) 
approximation of the amoeba. 



The diameter of inner complement components. We briefly discuss that the SOS- 
based certificates can also be used for more sophisticated questions rather than the pure 
membership problem. For this, we consider a family of polynomials / G C[Z\, . . . , Z n ] 
whose Newton polytope is a simplex and which have n + 2 monomials such that one of 
them is located in the interior of the simplex. Amoebas of polynomials of this class have 
at most one inner complement component (see [23] for a comprehensive investigation of 
that class). 

Let / := ^S) 1 hi ' Z a( - % \ and let bo denote the coefficient of the inner monomial. By 
the results in [25], for |6q| — > 00 the inner complement component appears at the image 



under the Log-map of a minimal point 5 of the function / :- 



/ 



5 is explicitly 



rg(fco)^ Q (°> 

computable, and \8\ is unique. This allows to certify that a complement component of 
the unlog amoeba has a certain diameter d under the scaling \Z\ i-> \Z\ 2 of the (unlog) 
amoeba basis space by solving the SDP corresponding to 

3 

J2^9i + H+1 = 
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Figure 5. The amoeba of / = X 2 Y + XY 2 + c • XY + 1 approximated 
with SOStools for c = 2 and c = -4. 

with Sl = ELM\ 2 - iZ^r-d'/A, s 2 = Eto (bi ■ Z^J e and s, = E^c? (* ■ Z*®)*", 
where gi G C[Z] (restricted to some total degree) and H is an SOS polynomial. 

Feasibility of the SDP certifies that there exists no point v G V(/) fl dBd/2(8) (where 
Bd/2{$) denotes the ball with radius d/2 centered in 5) in the rescaled amoeba basis space. 
Hence, the corresponding inner complement component of the unlog amoeba has at least 
a diameter d in that space. We have to investigate the rescaled basis space of the unlog 
amoeba in order to transform the generic condition (\5i\ — \Zi\) 2 — d 2 /4 on the basis space 
of Uf into a polynomial condition, which is given by s\ here. 

Note that this works not only for polynomials in the class under investigation, but for 
every polynomial as long as one knows, where a complement component appears. 

Example 5.3. Let / := 1 + Z\Z 2 + Z\Z\ + c- Z\Z 2 with a real parameter c. For this class, 
the inner complement component appears at the point (1,1) and thus under the Log- map 
at the origin of Log(M 2 ). The inner complement component exists for c > 1 and c < —3 
(cf. [25]). We compute a bound for the diameter of the inner complement component 
using the upper SDP for the intervals [—3,-9] and [1,7] with steplength 0.1. For any 
of these points we compute 14 SDPs in order to estimate the radius (based on binary 
search). In the rescaled amoeba basis space we obtain the bounds shown in Figure 

Observe that these bounds are lower bounds since feasibility of the SDP certifies mem- 
bership in the complement of the amoeba but infeasibility only certifies that no certificate 
with polynomials of degree at most k (i.e., 3 in our case) exists. 

This approach also yields lower bounds for the diameter of the inner complement com- 
ponent of the (log) amoeba. The image of the circle Xir=i(l^| 2 — \Zi\ 2 ) 2 — r<1 under rescaling 
and the Log-map, i.e. ^r=i(l°§ 1^1 — l^l) 2 — r2 5 contains the set of points 

[(\5 1 \-e r ,...AS n \),...A\Si\,---A^\-e r ),( K \5 1 \-e^,...,\S n \-e^)}. 
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Figure 6. Lower bound for the diameter of the inner complement component. 

By convexity of the complement components in Af, the simplex spanned by these points 
is contained in the inner complement component. Hence the double radius of the insphere 
of that simplex is a lower bound for the diameter of the inner complement component 
of A f . 

6. Outlook 

We have developed foundations and techniques for approximating amoebas and coamoe- 
bas based on the real Nullstellensatz and sums-of-squares techniques. While our focus was 
on developing the core principles of the computational methodology, some experimental 
results were provided to show the validity of the approach. Beyond the specific results 
we have presented, our approach can be seen as a first systematic treatment of amoebas 
from a (computational) real algebraic point of view and we think that this viewpoint will 
have more potential to offer. 

Major current challenges involve both computational issues (concerning the quality 
and efficiency of computations) as well theoretical questions on on the real algebraic 
viewpoint on amoebas. To name a specific open question in the latter respect, recall that 
in Theorem 14 A\ for a hypersurface amoeba Af we could deduce the order of a complement 
component from the special certificates we treated in that theorem. However, it is an open 
and important structural question how to deduce the order for a point w in the complement 
of Af given an arbitrary Nullstellensatz certificate for w. 

Acknowledgment. Thanks to Mihai Putinar for pointing out reference [Hj. Part of this 
work was done while the first author was visiting Universidad de Cantabria in Santander. 
He would like to express his gratitude for the hospitality. 
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